Universal topological phase of 2D stabilizer codes 
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Two topological phases are equivalent if they are connected by a local unitary transformation. In 
this sense, classifying topological phases amounts to classifying long-range entanglement patterns. 
We show that all 2D topological stabilizer codes are equivalent to several copies of one universal 
phase: Kitaev's topological code. Error correction benefits from the corresponding local mappings. 
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The theory of Ginzburg and Landau has had a tremen- 
dous success at classifying the different phases of mat- 
ter in terms of local order parameters and spontaneously 
broken symmetries. However, it fails to classify certain 
states of nature, such as the different fractional quantum 
Hall fluids which all have the same local symmetries. The 
Hamiltonian of these systems has a constant energy gap, 
and the ground state degeneracy depends on the topology 
of the space. Crucially, all ground states are locally iden- 
tical, which explains the failure of the Ginzburg-Landau 
paradigm. Instead, the classification of these systems re- 
quires the concept of topological order. 

Because topological order reflects the long-scale many- 
body correlations of the system, it cannot be modified 
locally. This robustness [T is indeed one of the many fea- 
tures that makes topologically ordered systems interest- 
ing for quantum information processing [2 . It also sug- 
gests a natural classification of topological phases: sys- 
tems that only differ by a local rearrangement of their de- 
grees of freedom belong to the same topological phase. In 
other words, the different phases are characterized only 
by their long-range entanglement patterns [3 . 

Another, more conventional, description of these 
phases is in terms of adiabatic connections. If two local 
and gapped Hamiltonians are connected by a family of 
local and gapped Hamiltonians, then it should be possi- 
ble to adiabatically interpolate between the two without 
encountering a phase transition. The two systems should 
therefore be in the same phase. This adiabatic evolution 
will generate a local unitary transformation [4] , so conse- 
quently the two systems will be in the same topological 
phase according to the definition adopted above. 

Quantum error-correcting codes [5 are intimately re- 
lated to topological order. To protect the information 
from local errors, information is encoded into the long- 
range entanglement of the system. A stabilizer code [6] 
is a special type of quantum code that can be defined 
as the degenerate ground state of a Hamiltonian on N 
qubits of the form 

H = -Y,^^ with [Sa,Sb]={) \/a,b (1) 

a 

where the stabilizer operators Sa are Hermitian elements 
of the Pauli group, i.e. they are constructed from tensor 



products of the three Pauli matrices cr^^, <j^, and cr^ and 
the identity operator /. Stabilizer codes are also frustra- 
tion free^ meaning that the Sa do not generate —1 under 
multiplication, so the ground states of H are +1 eigen- 
states of all stabilizers, i.e., S'al'0) = +1'^) for all a. The 
Sa form an Abelian group under multiplication, the sta- 
bilizer group S. When the qubits are embedded on a reg- 
ular lattice, the code — or its associated Hamiltonian — is 
said to be local if each operator Sa has support on a re- 
gion of constant size, independent of the system size. The 
support of an operator contains those qubits on which it 
acts nontrivially. 

In this Letter, we are interested in stabilizer codes that 
(i) are local and translationally invariant (LTI), and (ii) 
are topological, in the sense that no local operator can re- 
cover any encoded information — i.e., they have a macro- 
scopic minimum distance in terms of error correction. If 
we place our stabilizer in an infinite lattice, this can be 
formalized as follows. 

Definition 1 A topological stabilizer code (TSC) is a 
LTI stabilizer S such that Z{S) ex S. 

The symbol Z{S) denotes the centralizer of <S, the group 
of Pauli operators (with bounded support) that commute 
with all the elements of S. Our main result is that the 
topological phase of any 2D TSC is uniquely determined 
by its total quantum dimension or equivalent ly by its 
topological entanglement entropy ^topo = n\og2 [7, 8j. 
This follows from the existence of a local mapping io n/2 
copies of Kitaev's topological code (KTC) [10]. We 
also adapt the result to a class of subsystem stabilizer 
codes tlim. 

Many considerations motivate this line of research. 
Firstly, stabilizer codes provide simple models to study 
many-body quantum physics because they often admit 
exact solutions, and at the same time can exhibit com- 
plex phenomena such as topological order and anyonic 
excitations [9l[T0l[T3]. To our knowledge, this is the first 
example where the definition of topological order based 
on local equivalence [3 can be directly applied to a class 
of models in a rigorous manner. Secondly, in the con- 
text of error correction, the local equivalence to KTC 
enables us to directly extend a number of properties of 
this code to all 2D TSCs. For instance, thermal instabil- 
ity [Ml [15], code tradeoffs [16 , logical operator geometry 
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[17], and scale invariance [22] all become trivial corollar- 
ies of our mapping. In addition, our mapping provides a 
method to decode any 2D TSC code, while only a hand- 
ful of special cases previously had solutions [IHl UHl US • 
Thirdly, the local mapping can be used to change encod- 
ing during a quantum computation. Because the map- 
ping is local, this change will not propagate errors and 
is therefore fault-tolerant. This allows to put together 
the features of different codes — such as having transver- 
sal Clifford gates [13 , lower weight stabilizer generators 
[9l [TOl [20], etc. — and suggests a natural generalization 
of the notion of transversality for topological codes to 
include all local gates. 

Definitions — The notion of locality plays a crucial role 
here. For an operator X acting on the qubits of a 2D lat- 
tice, let us denote by |X| the range of X, defined as the 
size of the smallest square containing the support of X. 
With this definition, a Hamiltonian of the form Eq. ([T]) 
is local if there exists a constant w such that l^^l < w 
for all a. A translationally invariant unitary transfor- 
mation U is local if there exists a constant v such that 
\WXU\ < \X\ for ah operator X. Note that this 
definition is equivalent [21 to the requirement that U be 
decomposable into a system-size independent sequence 
of nearest neighbor unitary transformations. Lastly, we 
will say that two local stabilizer codes defined by Hamil- 
tonians H and Eq. ([T]), with stabilizer groups S and 
are locally equivalent if there exists a local unitary 
U and two trivial LTI stabilizer groups T and such 
that U {S ^ T)U^ = S' ^ T' . A trivial stabilizer group is 
generated only by single-qubit operators. Physically, U 
takes the ground state of H onto that of i^', and adds 
or removes extra qubits that are completely unentangled. 
The existence of renormalization group transformations 
that disentangle some qubit from topological codes [22] 
show the necessity of T and in this definition. 

Kitaev's topological code [9] [10] is defined on a 2D 
square lattice, with one qubit attached to each edge. For 
each lattice site 5, define an operator As = YI^^e 
where Eg denotes the set of edges incident to site s. Sim- 
ilarly, define for each lattice plaquette p (site of the dual 
lattice) an operator Bp = Yl^^^ cf% where Ep denotes 
the set of edges adjacent to plaquette p. The Hamilto- 
nian of the model is 

H = -Y,A,-Y,B,. (2) 

S p 

The excitations are anyons, gapped and topologically 
charged. Indeed, any set of excitations on the KTC can 
be reduced by local operations to one of four configura- 
tions: the vacuum (0) corresponding to no excitations, 
an electric charge (e) corresponding to a plaquette exci- 
tation Bp, a magnetic charge (m) corresponding to a site 
excitation As, and a composite excitation (/) containing 
both. These four sectors are the topological charges of 
the model. Excitations with different charges are char- 
acterized by different topological interactions or braiding 




FIG. 1: a) Regular square-octagon lattice for TCC. The dia- 
monds can be labeled A or B according to a chessboard pat- 
tern. There are two stabilizers Eq. ([s} associated to each pla- 
quette. b) Expanded square-octagon lattice for TSC. Starting 
with the lattice on the left, each vertex is expanded into a tri- 
angle. There is one gauge operator Eq. Q per edge, c) Zoom 
of a region of the extended lattice and rearrangement of the 
qubit into three stacks. 

statistics. According to the effect of exchanging two iden- 
tical charges, electric and magnetic particles are classi- 
fied as bosons, while the composite particle is a fermion. 
As for mutual statistics, they are all semionic because 
braiding any two distinct non-vacuum charges yields a 
— 1 phase. Finally, two charges can merge to form a new 
charge. The corresponding fusion rules are Abelian and 
such that m X e ^ f and a x a ^ for a = m, e, /. 
The notion of topological charge is of utmost relevance 
because local equivalence preserves the anyon model. 

Main result — We assume that 2D TSCs cannot give rise 
to chiral anyons [27 because the Hamiltonian terms Sa 
commute with each other [23]. Under this assumption, 
our main result is: 

Theorem 1 Every 2D TSC is locally equivalent to a fi- 
nite number of copies of KTC. 

By n copies of the code, we mean stacking n lattices on 
top of each other, each with the same Hamiltonian ([2|. 
This result implies that equivalence classes are labeled by 
the total quantum dimension of the code. The proof is 
rather technical and we do not present it here. It has two 
main steps. The first one shows that the excitations of a 
2D TSC have the same topological charges as a number 
of copies of KTC. Then, after a suitable renormaliza- 
tion, the hopping terms of the corresponding anyons are 
mapped to each other. 

We illustrate this for topological color codes (TCCs) 
[131 [20]. A TCC can be constructed on any 3-valent lat- 
tice with 3-colorable faces, but we take in particular the 
square-octagon regular lattice of Fig.[l]a). This lattice is 
particularly useful in terms of fault-tolerance [13] . Qubits 
are located at the vertices of the lattice, and there are two 
stabilizer operators per plaquette p 

= (8)^'' withaGK,a,}. (3) 

eeEp 

The excitations in this model carry 16 different topolog- 
ical charges that correspond exactly to the charges ob- 
tained from two copies of KTC. Guided by this charge 
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TOO-A KT01 KT02 TCC-A KTC1 KT02 TCC-B KTC1 KTC2 TCO-B KTC1 KTC2 



FIG. 2: Mapping between the Pauli operators of the square- 
octagon TCC an two copies of Kitaev's code KTCl, KTC2. 
The first (last) two columns are for the A (B) sub-lattice. 
Circles (stars) represent az (ctx) operators. For instance, the 
upper left diagram indicates that a ax located at the top 
of a diamond of the A sub-lattice gets mapped to a ax on 
KTCl and two az on KTC2. All commutation relations are 
preserved by this mapping, so it is unitary and obviously local. 



identification, we obtain the mapping shown at Fig. |2] It 
can be directly verified that it maps stabilizer generators 
of TCC to stabilizer generators of two KTC, in this case 
with no need to add trivial stabilizers. 

Subsystem codes — Subsystem stabilizer codes form a 
more general class of stabilizer codes E] . They can 
be defined as a pair (5, Q), where Q is an arbitrary Pauli 
subgroup and S a stabilizer such that S oc Z{Q) f] Q. 
Encoding is not done on the whole subspace defined by 
5, but rather on the subsystem where the action of Q is 
trivial. This way, errors caused by operators in Q do not 
affect encoded states. Because of this, elements of Q are 
called gauge operators. 

We say that a subsystem code (S^Q) is LTI if Q admits 
a LTI set of generators Q^. Note that some local subsys- 
tem codes admit no local stabilizer generators, e.g. [12 . 
Unlike them, a topological subsystem code should have 
a stabilizer with a local description. In addition, local 
operators should not recover any encoded information. 
Since we do not care about the effect of gauge operators, 
this can be formalized as follows in an infinite lattice: 

Definition 2 A topological stabilizer subsystem code 
(TSSC) is a LTI subsystem stabilizer code {S^Q) such 
that Z{S) (X Q. 

There is a general strategy to understand TSSCs in 
terms of TSCs. Namely, to find a TSC S' that lies in 
between the stabilizer group and the gauge group of the 
subsystem code, i.e. S d S' d Q. We can then map S' 
invoking Theorem 1, which shows that S is locally equiv- 
alent to a subset of the generators of several copies of 
KTCs. Not all topological charges of these KTCs carry 
information: some are associated to gauge operators. 
These "gauge charges" do not topologically interact with 
charges that describe errors in that we call "proper 
charges". Stabilizer of S detect the presence of proper 
charges, ignoring any gauge charge. Unlike for TSCs, 
proper charges can give rise to a chiral anyon model as 
exemplified below. 



Let us illustrate this strategy with an important family 
of 2D subsystem codes [20 called topological subsystem 
color codes (TSCCs). Given the lattice of a TCC, we 
can inflate each vertex into a triangle as in Fig. [l]b). 
Qubits are located on the vertices of this inflated lattice, 
and there is one gauge group generator associated to each 
pair of sites j connected by an edge 
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with a = cFx-,cFy^ OT cFz for a dashed, dotted, or solid edges 
respectively. This code admits a set of local stabilizer 
generators, some of which involve a relatively large num- 
ber of qubits (up to 24). Errors are described by three 
nontrivial topological charges, all fermions (/i,/2,/3), 
making it a chiral anyon model. All the mutual statis- 
tics are semionic and the fusion rules are / x / ^ and 
fi X fj fk with i, j, /c all different. These charges can be 
obtained as a subset of two copies of KTC. For instance, 
one can identify fi ^ [m, /], /2 ^ [e, /], and /s ^ [/, 0] 
where the notation [a, b] stands for a composite particle 
of charge a on the first KTC and charge b on the second. 
This suggests that we should be able to find a TSC 
with S C S' C Q and locally equivalent to n copies 
of KTC with n > 2. We will present two different ways 
of obtaining S' that are geometry independent (i.e., not 
restricted to the square-octagon lattice). 

In the first construction, S' is the stabilizer of three 
TCC on the corresponding non-inflated lattice. Indeed, 
all we need to do is to rearrange the qubits. The three 
qubits located at the vertices of each triangle inherit the 
color label of the neighboring plaquette. We construct 
a stack of three TCC lattices — one per color — each one 
containing the qubits of that color, see Fig. [l]c). It can 
be easily verified that this maps the generators of 5 to a 
subset of the generators of the three TCCs. We obtain S' 
by including the other generators of these TCCs. In the 
second construction, we consider the stabilizer group Sz 
generated by the gauge operators of the form a^cr^ (solid 
edges), which clearly is a subgroup of Z{S). It follows 
from the results in [24 that S' := SSz is a TSC with 
the same topological charges as a TCC. These two con- 
structions illustrate that the quantum dimension of the 
intermediate code S' is not uniquely determined, since in 
the first case we have Hi = 2^ and in the second = 2^. 

Because 5 is a strict subset of S' ^ the mappings de- 
scribe above do not take the system to the ground state 
of the resulting KTCs; the stabilizer added to S to arrive 
at will generally contain excitations. These can be 
eliminated by local transformations, fusing particle pairs 
into the vacuum. Moreover, because these excitations 
correspond to gauge charges, this local transformation 
does not change the encoded information. 

Decoding — When the system is prepared in the ground 
state of the Hamiltonian Eq. ([T]), all stabilizers have value 
+1. But in the presence of errors, this will not be the 
case in general. The problem of decoding a quantum code 
consists in identifying the most likely recovery to restore 
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FIG. 3: Decoding failure probability as a function of the error 
probability of each qubit for the octagon-square TCC (left) 
and TSCC (right), based on the algorithm of [H]. The dif- 
ferent curves illustrate lattices of different linear size /: below 
a threshold probability (dotted lines), the decoding failure 
probability decreases with the lattice size, leading to a per- 
fect recovery in the thermodynamic limit. 



the encoded state from partial information coming from 
the measurement of the stabilizer operators, whose ±1 
outcomes are called error syndrome. Not all codes can 
be decoded efficiently, but fast approximate algorithms 
have been devised for KTC. The one presented in [10] 
uses Edmonds' minimum matching algorithm [25^ to find 
the shortest path that recombines all electric particles in 
pairs and independently all magnetic particles in pairs. 
For N qbits, it runs in time N^. The algorithm pro- 
posed in [18 uses renormalization group approximations 
to find the homological class of errors with the highest 
probability. It runs in time logA^. An efficient decoder 
was also devised for a TSSC on a particular lattice [19], 



but in general each new code requires a tailored decoding 
technique. 

Our techniques can be used to decode any 2D TSCs 
or TSCCs. The idea is to choose a set of "elementary" 
charges — either fermionic or bosonic — that generate all 
the topological charges in the code. Then, the decod- 
ing problem can be mapped to that of KTC described 
above for each elementary charge. We have used this 
technique, combined to the decoding algorithm of [18], 
for the TCC on the square-octagon lattice of Fig. [T] on a 
bit-flip channel and found an error threshold of roughly 
8.7%, in good agreement with the Monte Carlo estimate 
of 10.9% [26^ for ideal error correction. We have also 
used this technique for the TSCC on the square-octagon 
lattice of Fig. [T] on a depolarizing channel and found an 
error threshold of roughly 1.95%, in good agreement with 
the estimate of 2% [19] for a closely related code in a five- 
square lattice. 

Conclusion — We have demonstrated that 2D topologi- 
cal stabilizer codes all belong to one universal topologi- 
cal phase by constructing an explicit local mapping onto 
multiple copies of Kitaev's topological code. This results 
also carries to a certain class of 2D subsystem codes, 
and in particular to all topological subsystem color codes. 
These local maps enable us to extend many properties of 
Kitaev's code to all 2D codes, and in particular directly 
yield efficient decoding algorithms for error correction. It 
could also have important implications for fault-tolerant 
quantum computation. 
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